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SUMMARY 


The  results  of  a numerical  and  experimental  investigation  of 


i 


friction  and  heat  transfer  parameters  for  turbulent  flow  of  helium-argon 
mixtures  in  smooth,  electrically  heated,  vertical  circular  tubes  are 
presented.  Mixtures  with  molecular  weights  between  15.3  and  29.7  are 
used;  these  values  correspond  to  molecular  Prandtl  numbers  between  0.42 
and  0.49.  Inlet  Reynolds  numbers  range  from  31,200  to  102,000,  maximum 
wall  tenperatures  from  392  to  828°K,  maximum  wall-to-bulk  temperature 
ratios  to  1.82,  and  pressures  from  4.7  to  9.7  atmospheres.  Existing 
experimental  correlations,  developed  using  gases  with  Prandtl  numbers  of 
the  order  of  0.7,  are  found  to  overpredict  the  observed  Nusselt  numbers. 

By  comparison  of  numerical  calculations  and  measured  constant  property 
Nusselt  numbers , turbulent  Prandtl  numbers  are  determined  in  the  wall 
region.  The  validity  of  using  these  deduced  turbulent  Prandtl  numbers  is 
examined  for  conditions  where  the  properties  vary  significantly. 

Heat  trains fer  and  wall  friction  parameters  have  been  obtained 
numerically  to  demonstrate  the  effects  of  mixture  composition  and  gas 
property  variation  for  heating  or  cooling  in  regenerative  heat  exchangers . 
The  situation  is  modelled  by  laminar  flow  through  short  ducts  with  constant 
wall  heat  flux.  For  design  predictions  accounting  for  the  effect  of  prop- 
erty variation,  it  is  found  that  the  property  ratio  method  is  better  than 
the  film  temperature  method  for  heat  transfer  while  the  latter  method  is 
preferable  for  apparent  wall  friction  - with  the  proviso  that  the  present 
definitions  of  the  non-dimensional  parameters  be  employed. 

iii 


The  present  status  of  experiments  using  alternate  gas  mixtures 
and  examining  the  effects  of  heating  rate  on  heat  transfer  from  rough 

surfaces,  such  as  ceramics  or  artificially  roughened  superalloys,  is  pre 
sented  briefly. 


NOMENCLATURE 


a exponent  used  to  account  for  temperature  variation  of  viscosity; 

A cross  sectional  area  of  tube; 

cs 

b exponent  used  to  account  for  temperature  variation  of  thermal 

conducti vi ty ; 

c velocity  of  sound; 

c^  specific  heat  at  constant  pressure; 

D inside  diameter; 

g gravitational  constant; 

gc  dimensional  conversion  factor; 

G average  mass  flux,  m/A  ; 

cs 

h heat  transfer  coefficient; 

i specific  enthalpy; 

k thermal  conductivity; 

l mixing  length ; 

m mass  flow  rate; 

M molal  mass ; 

p pressure; 

q"  heat  flux; 

r radius ; 

R gas  constant  for  a particular  gas; 

R universal  gas  constant; 

T temperature ; 

u velocity  in  axial  direction; 

v velocity  in  radial  direction; 


v 


driu/tm«-tiir‘it  "i  m 


x axial  distance  from  start  of  heating; 

y radial  distance  from  wall. 

Greek  Symbols 

£/<  force  constant  in  Lennard- Jones  potential; 

eddy  diffusivity  for  heat; 

£m  eddy  diffusivity  for  momentum; 

Y ratio  of  specific  heats,  c^/c^; 

< empirical  constant  in  van  Driest  mixing  length  model,  0.4 

y absolute  viscosity; 

V kinematic  viscosity; 

p density; 

O'  force  constant  in  Lennard- Jones  potential; 

t shear  stress. 


Non-dimensional  Parameters 


Gr 

L+ 

L* 

Nu 

P 

Pr 

+ 

q 

Q+ 

Re 


friction  factor,  2gcPT ^/G  ; 

4 2 

Grashof  number  based  on  wall  heat  flux,  gD  q"  /(v  yc  T) . j 

w pi 

length  for  velocity  boundary  layer,  4L/(D^Re) 

length  for  thermal  boundary  layer,  4L/(D  RePr) 

h 

Nusselt  number,  hD/k; 

2 

pressure  drop,  p^g^p^pJ/G  ; 

Prandtl  number,  c yA;  Pr..,  turbulent  Prandtl  number; 

p t 

heat  flux  parameter,  q£/(Gc^, /TJ  ; 
laminar  heat  flux  parameter,  q^D^/(k^T^) ; 

Reynolds  number,  GD/y; 


y+  wall  distance  parameter,  y (gcT^p)  1//2/V; 

y+  empirical  constant  in  van  Driest  mixing  length  model,  26. 

Subscripts 

av  lengthwise  average; 

b evaluated  at  bulk  temperature ; 

cp  constant  property  condition; 

DB  Dittus-Boelter; 

DKM  Drew,  Koo  and  McAdams; 

eff  effective; 

i inlet;  an  index; 

Max  maximum ; 

ref  reference; 

t turbulent ; 

VD  van  Driest; 

w wall. 
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I.  INTRODUCTION 


1.1.  Background 


The  closed  Brayton  cycle,  or  gas  turbine  cycle,  has  been  suggested 
as  an  efficient,  compact,  versatile  system  for  power  plant  and  propulsion 
applications  [Bammert,  Rurik  and  Griepentrog,  1974;  Mock,  1970].  With 
thermodynamic  cycle  studies,  Bammert  and  Klein  [1974]  showed  that  con- 
siderable savings  in  the  total  costs  of  the  turbomachines  and  heat  ex- 
changers can  be  achieved  by  mixing  helium  with  a heavier  gas.  When  the 
heavier  gas  is  another  noble  gas , approximate  calculations  indicate 
the  possibility  of  significant  improvements  in  heat  transfer  performance 
compared  to  pure  gases  at  the  same  conditions  [Vanco,  1965].  However,  the 
gaseous  data  for  heat  exchangers  and  analyses  'or  idealized,  but  related, 
conditions  have  concentrated  on  pure  gases,  primarily  air,  with  constant 
transport  properties.  Whether  such  results  can  be  applied  for  mixtures 
of  inert  gases  with  temperature-dependent  properties  is  a basic  question. 

The  closed  Brayton  cycle  offers  advantages  of  reliability,  high 
thermal  efficiency,  reduced  atmospheric  pollution,  minimal  erosion  or 
blockage  of  components,  rapid  startup  and  high  part-load  efficiency  to  a 
wide  range  of  uses.  It  has  been  proposed  for  application  to  surface  and 
undersea  propulsion  of  naval  ships  in  conjunction  with  gas-cooled  nuclear 
reactors  and  fossil  fueled  combustors,  and  for  auxiliary  electrical  systems, 
for  solar  power,  power  systems  for  satellites  and  space  stations,  bus  and 
rail  power  plants  and  refrigeration/heating  systems.  Cycle  efficiency  is 


expected  to  improve  as  technology  advances.  By  employing  noble  gas  mixtures 
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such  as  helium  and  xenon  as  the  working  fluid,  heat  transfer  performance 
can  be  improved  and  system  volume  and  weight  can  be  reduced  compared  to 
common  gases  at  the  same  conditions.  The  understanding  and  effective  use 
of  these  methods  of  increasing  heat  transfer  performance  are  important  to 
the  design  of  lightweight  propulsion  plants. 

1.2.  Conditions  In  Naval  Brayton  Systems 

The  non-dimensional  governing  parameters  have  been  estimated  for 
operating  conditions  in  the  Brayton  test  unit  at  the  Naval  Ship  Research 
and  Development  Center.  With  argon  as  the  working  fluid,  the  following 
values  are  expected  (approximately) . 


Unit 

Reynolds 

Number 

Grashof 

Number 

(buoyancy) 

Mach  Number 
(Compress- 
ibility) 

Wall-to-Bulk 
Temp . Ratio 

Heater  Tubes 

50,000 

560 

0.3 

* 

1.15 

Recuperator 

Low  Pressure 

1,800 

110 

0.06 

0.9 

High  Pressure 

2,700 

160 

0.05 

1.1 

Cooler 

2,300 

2,000 

0.04 

0.85 

Thus,  buoyancy  effects  are  expected  to  be  unimportant.  The  heater  tubes 
would  have  fully  turbulent  flow  and,  in  further  modification  of  the  design, 
compressibility  effects  could  become  important.  The  Reynolds  number  range 
in  heat  exchangers  implies  either  laminar  or  transitional  flow.  Laminar 
flow  might  be  predicted  well  with  a three-dimensional  numerical  analysis 
but  the  transition  region  would  require  experiments.  The  temperature  ratios 
indicate  modest  property  variation  in  present  designs  but  with  a graphite 
heater  or  ceramic  components  extreme  property  variation  could  become  involved. 
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1 . 3.  Proposed  Work 


The  initial  proposal  and  subsequent  proposals  by  the  University 
of  Arizona  [McEligot,  1974a, b,  1975a, b]  detailed  the  applications,  advant- 
ages of  inert  gas  mixtures,  convective  problems,  prediction  of  heat  trans- 
fer and  pressure  drop,  scientific  needs,  previous  work  at  the  University 
of  Arizona,  proposed  experimental  and  analytical  studies  and  related  work, 
facilities  and  apparatus  available,  personnel,  and  nomenclature.  The 
general  goal  is  to  develop  the  ability  to  predict  accurately  temperature 
and  velocity  distributions , pressure  drop  and  wall  heat  fluxes  for  flows 
in  components.  Scientific  needs  related  to  the  proposed  work  were: 

1.  Determine  whether  existing  scaling  laws  are  adequate  for 
mixtures  of  inert  gases. 

2.  Extend  numerical  analyses  to  include  inert  gas  mixtures  in 
the  range  0.2  <!  Pr  < (5.7. 

3.  Test  the  analyses  experimentally. 

4.  Obtain  definitive  data  on  the  effect  of  molecular  Prandtl 
number  on  turbulence  models. 

5.  Obtain  greater  understanding  of  turbulent  shear  flow  through 
the  proposed  studies. 

6.  Determine  the  effects  of  element  geometry,  Prandtl  number  and 
heating  rate  on  convective  heat  transfer  from  rough  surfaces . 

Accordingly , objectives  were  summarized  as : 

1.  For  mixtures  of  inert  gases  with  temperature-dependent  trans- 
port properties,  determine  the  effects  of  mixture  composition 
and  heating  rate  on: 

a.  convective  heat  transfer  parameters,  such  as  the  Nusselt 
number  and/or  the  Stanton  number, 

b.  pressure  drop  or  wall  friction,  in  terms  of  the  friction 
factor, 
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c.  convective  heat  transfer  and  fluid  friction  for  flow  over 
artificially  roughened  surfaces  in  terms  of  the  roughness 
functions . 

2.  Extend  theoretical  predictions  into  range  0.2  < Pr  < 0.7  and 
improve  understanding  of  turbulent  heat  transfer  with  property 
variations  for  both  smooth  and  artificially  roughened  surfaces. 

Tasks  suggested  include: 

1.  Numerical  predictions  for  inert  gas  mixtures  in  heated 
circular  tubes  as  in  some  proposed  heaters  for  ship  propulsion. 

2.  Preliminary  measurements  of  pure  gases  and  mixtures  in  open- 
loop  apparatus;  comparison  to  predictions  (first  and  second 
years) . 

3.  Design  and  modifications  of  apparatus  to  provide  closed  loop 
operation  with  the  more  expensive  mixtures;  measurements  with 
pure  gases  and  mixtures;  comparison  to  predictions. 

4.  Design  and  fabrication  of  test  sections  to  obtain  mean  velocity 
and  temperature  profiles  for  flow  over  artificially  roughened 
surfaces  with  heat  addition. 


2.  HEATED  TURBULENT  FLOW  OF  HELIUM-ARGON 
MIXTURES  IN  TUBES 


2.1 . Introduction 

In  preliminary  design  studies  for  application  of  inert  gas  mixtures 
to  closed  gas  turbine  systems,  Vanco  ^19653  estimated  relative  heat  transfer 
coefficients  for  turbulent  flow  in  ducts  using  a version  of  the  Colburn 
analogy  CKreith,  19583  which  may  be  written  as: 

Nu  = 0.023  Re°*8  Pr1/3  (1> 

This  relationship  is  supposedly  valid  for  the  range  0.5  < Pr  < 100,  but  is 
primarily  based  on  experiments  for  fluids  with  Pr  ~ 0.7  or  greater. 


For 


gases,  with  Pr  - 0.7,  McAdams  Cl9543  recommended  the  Dittus-Boelter 
correlation  with  the  coefficient  modified 

Nuqb  = 0.021  Re0'8  Pr°'4  (2) 

instead.  The  difference  is  about  ten  percent  for  most  gases  (with 
Pr  ~ 0.7 J.  As  3hown  later  the  Prandtl  number  for  a mixture  of  helium 
and  argon  can  be  of  the  order  of  0.4;  no  data  are  available  for  this 
range  so  it  is  important  to  test  the  validity  of  these  correlations.  This 
testing  is  one  of  the  purposes  of  the  present  experiment. 

For  predictions  of  heat  transfer  to  turbulent  flows  many  engineer- 
ing analyses  employ  the  turbulent  Prandtl  number,  Pr  = e /e,  , to  estimate 

t m n 

effective  thermal  conductivity  from  correlations  for  momentum  transfer. 

As  shown  in  surveys  by  Blom  Ql970],  by  Quarmby  and  Quirk-*  Cl972]  and  by 
A.  J.  Reynolds  Cl975],  there  still  exists  considerable  uncertainty  con- 
cerning the  proper  variation  of  Pr  with  molecular  Prandtl  number.  At 
Pr  Z 0.7,  some  idealized  analyses  suggest  that  Pr  > 1 while  others  yield 
Pr  <1.  A second  purpose  of  the  present  work  is  to  determine  whether  Pr 
differs  significantly  from  unity  at  Pr  ~ 0.4  - 0.5  in  the  wall  region, 
which  dominates  the  convective  heat  transfer  behavior. 

The  next  section  summarizes  pertinent  knowledge  of  the  transport 
properties  of  mixtures  of  helium  and  argon.  Section  3 briefly  describes 
the  numerical  analysis  to  be  employed  and,  for  fully  established  conditions, 
presents  the  predicted  Nusselt  numbers,  which  disagree  with  correlations 
(1)  and  (2)  when  Reynolds  analogy  (Pr  = 1)  is  employed.  The  experiment  is 
then  described.  In  Section  5 the  data  are  extrapolated  to  the  constant 
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property  idealization  to  test  the  correlations  and  then  the  diagnostic 

technique  of  McEligot,  Pickett  and  Taylor  Cl976]  is  applied  to  estimate 

Pr  . The  following  section  extends  the  data  to  heating  rates  where  the 
t / w 

temperature  dependence  of  the  transport  properties  becomes  significant 
and  demonstrates  that  numerical  predictions  based  on  the  value  of  Prfc 
deduced  from  measurements  at  low  heating  rates  are  still  adequate. 

2.2.  Transport  Properties  of  Helium-Argon  Mixtures 

The  properties  needed  for  this  study  were  compressibility,  viscosity 
thermal  conductivity,  specific  heat,  enthalpy,  speed  of  sound,  and  the  gas 
constant.  The  properties  of  air  have  been  studied  extensively;  for  the 
comparison  data  the  NBS  tables  of  Hilsenrath  et  al.Cl955]  were  used  in  this 
investigation.  The  properties  of  helium  and  helium- argon  mixtures  were 

•a 

calculated  theoretically.  For  all  gases  the  viscosity  and  thermal  con- 
ductivity were  assumed  to  be  independent  of  pressure. 

The  helium  and  helium-argon  mixtures  were  assumed  to  be  ideal 
gases,  thus  making  the  compressibility  equal  to  unity.  This  assumption  is 
reasonable  for  the  range  of  pressures  (100  - 970  kPa)  and  temperatures 
(290  - 830°K)  used  in  this  experiment.  Since  helium  and  argon  are  mona- 
tomic and  the  temperature  range  in  this  study  was  not  too  large,  the 
relation  ^Reynolds , 1968^ 


c = (5/2)  R = (5/2)R/M 
P 


was  used  to  calculate  the  specific  heat. 
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Using  the  ideal  gas  and  constant  specific  heat  assumptions,  one 
may  derive  simple  equations  for  the  enthalpy  and  speed  of  sound  ^Reynolds 
and  Perkins,  19683 


1 = c 


(T 


T ) 
ref 


and 


(4) 


c 


/ (5/3) RT 


The  Lennard-Jones  (6-12)  potential  can  be  employed  in  the  Chapman- 
Enskog  kinetic  theory  to  predict  thermal  conductivity,  viscosity  and 
Prandtl  number  of  binary  mixtures  of  inert  gases  QHirschf elder , Curtiss 
and  Bird,  19643.  There  has  been  considerable  experimental  study  of  the 
pure  gases  but.,  few  data  exist  on  the  mixtures. 

With  force  constants,  £/K  and  o,  suggested  by  Hirschf elder , 

Curtiss  and  Bird  the  predicted  viscosity  for  helium  falls  about  eight 
percent  below  the  data  of  Dawe  and  Smith  El97o3  and  Kalelkar  and  Kestin 
Cl9703  at  temperatures  around  900°C.  Likewise,  the  predicted  thermal  con- 
ductivity is  about  nine  percent  lower  than  the  measurements  of  Saxena  and 
Saxena  Cl9683  up  to  1100°C.  Similar  discrepancies  exist  for  argon.  Using 
force  constants  suggested  by  DiPippo  and  Kestin  Cl9693  instead  leads  to 
essential  agreement  with  the  values  recommended  by  the  Thermophysical 
Properties  Research  Center  [ Touloukian  and  Ho,  19703- 

The  mixture  properties  are  shown  in  Figure  1.  The  solid  curves 
were  calculated  with  the  force  constants  recommended  by  DiPippo  and  Kestin: 
a = 2.158A  and  £/<  = 86.2  K for  helium  and  0 = 3.292A  and  £/<  = 152.75°K 
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for  argon.  The  dashed  curves  compare  values  calculated  from  the  recom- 
mendations of  Hirschf elder , Curtiss  and  Bird. 

The  viscosity  of  the  mixture  varies  only  slightly  with  molecular 
weight  (composition) . From  pure  helium  to  pure  argon  the  variation  is 
only  about  fifteen  percent;  for  mole  fractions  of  argon  greater  than 
0.25  the  viscosity  is  almost  constant.  As  with  pure  gases  the  viscosity 

increases  with  temperature.  The  predicted  values  agree  well  with  the 

■ 

data  of  Tanzler  CTouloukian  and  Ho,  1970H  at  75  atmospheres  and  of  Kalelkar 
and  Kestin  Cl970]. 

The  mixture  thermal  conductivity  decreases  by  a factor  of  five  or 
more  as  the  molecular  weight  increases  from  pure  helium  to  pure  argon  and 
it  increases  with  temperature.  Agreement  with  the  data  of  Gambhir  and 
Saxena  Cl966l  and  of  Cheung,  Bromley  and  Wilke  Cl962]  is  close.  The 
measurements  of  Peterson,  Hahn  and  Commings  [l97l]  at  50  atmospheres  show 
that  neglect  of  pressure  variation  is  justified.  It  appears  that  the  higher 
temperature  data  of  Mason  and  von  ubisch  Cl960^  refute  the  model;  however, 
in  a critical  review  Gandhi  and  Saxena  Cl958^  have  observed  that  the  mea- 
surements of  von  Ubisch  appear  to  be  systematically  higher  than  others 
they  reviewed. 


In  a companion  study  for  laminar  flow,  McEligot,  Taylor  and  D irst 
Cl977j  demonstrated  that  the  temperature  dependence  of  the  viscosity  and 
thermal  conductivity  could  be  approximated  as 


The  exponent  a ranges  from  0.7  to  0.8  and  b falls  between  about  0.7  and 


0.75. 


As  a consequence  of  the  variation  of  thermal  conductivity  and 


specific  heat  versus  molecular  weight  the  Prandtl  number  decreases  to  a 
minimum  of  about  0.42  at  $ ~ 16  from  about  2/3  for  the  pure  gases.  It 
is  about  0.48  at  the  molecular  weight  of  air.  The  temperature  dependence 
is  almost  negligible  since  the  power  law  exponents  for  y and  k differ  so 
little. 

In  the  present  study  the  calculated  properties  were  applied  in 
tabular  form  in  the  data  reduction  programs.  For  numerical  analyses, 
correlations  of  the  calculated  values — such  as  equations  (5) — were  used 
to  reduce  computation  time. 


2.3.  Numerical  Analysis 


In  this  paper,  analysis  is  by  the  numerical  method  of  Bankston 
and  McEligot  Cl970]  utilizing  finite  control  volume  approximations  to 
solve  the  governing  equations. 

Continuity: 


3 3u  , 1 3 ao  _ 
3x  r 3r 

x-momentum : 
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Energy : 
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Integral  continuity: 


f (3u  rdr  = 1/2 
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Idealizations  implied  by  these  forms  are  (a)  the  axisymmetric  boundary 
layer  approximations,  (b)  steady  flow  at  low  velocities,  (c)  constant 
mixture  concentration,  and  (d)  negligible  axial  conduction.  The  circumflex 
C'')  represents  non-dimensionalization  with  respect  to  the  value  of  the 
quantity  of  the  entrance. 

Gas  properties  may  be  idealized  as 


a a 
P = p/T; 


Ta;  ft  = A- 


/v  £d 
C = T 
P 


Initial  conditions  are  specified  and  boundary  conditions  are  the  no-slip, 
impermeable  wall  with  the  observed  wall  heat  flux  variation  specified. 

The  effective  viscosity  is  predicted  with  a combination  of  the 
van  Driest  mixing  length  [1956]  and  Reichardt  middle  law  Ll95lJ, 


*Vd  = Ky  ^ " exP(_  y+/yA+) } 


e -e  _ . (2  - *-)  • Ci  + 2 (— ) 3/6 

m vD  r r 

w w 
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The  wall  temperature  is  used  to  evaluate  properties  in  y+.  With  < - 0.4 
and  y^+  = 26,  this  representation  yields  adiabatic  friction  factors  with- 
in about  one  percent  of  the  Drew  , Koo  and  McAdams  correlation  Cl9323  over 

4 5 

the  range  3 x 10  < Re  < 3 x 10  . 

The  program  uses  implicit  algebraic  equations  to  represent  the 
governing  equations.  These  equations  are  iterated  at  each  axial  step  to 
treat  their  coupling  and  the  non-linear  terms.  Mesh  spacing  is  chosen  to 
provide  Nusselt  numbers  and  friction  factors  within  about  two  percent  of 
their  converged  values.  The  node  nearest  the  wall  falls  at  y+  ~ 0.5  or  less. 

For  predictions  under  the  constant  property  idealizations  the  ex- 
ponents a,  b and  d are  set  to  zero  and  3 is  taken  as  unity.  Then  if  the 
inlet  condition  is  a fully  developed  flow,  only  the  energy  equation  is  solved 

Heat  transfer  results  for  fully  established  conditions  with  constant 
properties  are  shown  in  Figure  2.  Reynolds  analogy,  Prfc  = 1,  has  been  as- 
sumed here.  Normalization  by  the  Dittus-Boelter  correlation  (2)  shows  that 
such  correlations  may  be  expected  to  fail  as  the  Prandtl  number  is  reduced 
and  the  Reynolds  number  is  increased.  The  two  correlations  would  over- 
predict the  Nusselt  number  with  the  Colburn  analogy,  used  by  Vanco  Cl965l, 
being  about  fifteen  percent  worse  than  the  Dittus-Boelter  correlation  for 
helium-argon  mixtures.  Whether  Reynolds  analogy  and  these  numerical  pre- 
dictions are  adequate  must  be  determined  by  experiment. 

2.4.  Experiment 

The  experimental  apparatus,  arrangement,  and  procedure  were  similar 
to  those  used  by  Campbell  and  Perkins  [^1968  j.  Instead  of  a square  duct, 
a circular  tube  of  Hastelloy— X was  used  as  the  test  section.  This 
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ma 


Colburn  Pr  = 0.45 


vertical  tube  had  an  inside  diameter  of  0.312  cm  (1/8  in)  and  a wall 
thickness  of  0.056  cm  (0.022  in) . The  test  section  consisted  of  a heated 
section  98  diameters  in  length  preceded  by  an  unheated  section  92  diameters 
in  length.  The  onheated  section  ensured  that  the  velocity  profile  approached 
fully  developed  flow  prior  to  heating.  Electrical  resistance  heating  of 
the  test  section  yielded  an  axial  heat  flux  distribution  which  exponentially 
approached  a constant  value  within  three  diameters  and  then  remained  con- 
stant within  a few  percent.  Two  pressure  tapes  were  used.  One  was  located 
in  the  lower  electrode  and  the  other  eight  diameters  below  the  upper 
electrode. 

Sixteen  premium  grade  Chromel-Alumel  thermocouples,  0.013  cm. 

(0.005  in)  diameter,  were  spot  welded  to  the  heated  section  of  the  tube 
using  the  parallel  junction  suggested  by  Moen  Cl960^.  Thermocouple  con- 
duction error  was  calculated  from  the  heat  loss  calibration  data  and  a 
relation  developed  by  Hess  Ql965];  the  thermocouple  conductance  was  esti- 
mated from  the  emissivity  of  the  bare  wires  and  a natural  convection  cor- 
relation for  small  Rayleigh  numbers.  The  correction  is  of  the  order  of 
one  percent  of  the  difference  between  the  tube  temperature  and  the  en- 

_ 4 

vironmental  temperature;  at  Re  ~ 3 x 10  , the  effects  are  approximately 
I5!  to  2 percent  on  the  Nusselt  number  and  3 to  4 percent  on  the  deduced 
turbulent  Prandtl  number. 

With  the  exception  of  air,  the  gas  was  supplied  in  commercial  gas 
cylinders.  The  manufacturer,  Matheson  Gas  Products  Division,  provided  the 
mixtures  to  our  specifications  in  lots  of  four  bottles.  For  the  first  two 
batches  concentration  was  determined  by  gas  chromotography ; samples  of  the 


same  mixture  varied  by  two  percent  by  volume  so  adiabatic  friction  factor 
measurements  were  used  to  confirm  the  concentration  selected.  For  the 
last  two  batches  Matheson  determined  the  concentration  with  a thermal  con- 
ductivity column  with  sample  variation  to  within  0.8  and  0.4  percent  (how- 
ever, the  experimental  uncertainty  of  the  technique  is  estimated  as  two 
percent) ; these  concentrations  were  also  confirmed  by  the  adiabatic 
friction  data. 

To  measure  the  higher  flow  rates,  our  positive  displacement  meter 
was  replaced  by  a Meriam  laminar  flow  element.  The  latter  was  calibrated 
to  measure  the  flow  rate  within  + 1.5  percent.  Heise  gages,  inclined 
water  manometers,  and  vertical  mercury  or  water  manometers  were  used  to 
measure  static  pressure  and  pressure  drop. 

Table  1 summarizes  the  range  of  variables  covered  in  this  investi- 
gation. A more  detailed  discussion  of  the  experiment  and  tabulations  of 
the  data  are  available  in  a report  by  Pickett  ^1976^. 

The  experimental  uncertainties  were  estimated  by  the  method  of 
Kline  and  McClintock  Cl953].  Typical  values  for  the  Nusselt  number  are 
13  percent  at  x/D  ~ 1.2  decreasing  to  4 percent  at  x/D  ~ 25  for  low  heat- 


ing rates  and  slightly  less  for  the  higher  heating  rates . The  dominant 
uncertainty  is  in  the  bulk  stagnation  temperature  which  depends  on  the 


mass  flow  rate,  electrical  power  and  the  heat  loss  calibration. 

The  test  section  was  a bare  tube  surrounded  by  a draft  shield  so 
heat  loss  was  by  radiation  and  natural  convection.  The  heat  loss  was 
calibrated  as  a function  of  axial  position  and  temperature  from  measurements 
without  internal  flow.  During  runs  with  flow  the  fraction  of  power  dissi- 
pated to  the  environment,  instead  of  to  the  internal  gas  flow,  increases 
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Table  1.  Range  of  Variables  in  the  Present  Experiment 


Air 

He  1 i urn 

He  1 i um-Argon 

Experimental  Runs 

25 

4 

28 

Molecular  Weight 

23.97 

4.003 

15.3  - 29.7 

Inlet  Bulk  Reynolds  Number 

32,900-100,000 

30,200 

31,200-102,000 

Exit  Bulk  Reynolds  Number 

19,900-  89,000 

13,400-26,600 

17,000-  68,000 

Inlet  Bulk  Prandtl  Number 

0.719 

0.667 

0.419-0.486 

Exit  Bulk  Prandtl  Number 

0 . 602 

0.667 

0.426-0.495 

Maximum  T /T. 

w b 

1 .90 

1 .75  - 

1 .82 

Maximum  T (°K) 
w 

317 

739 

328 

Maximum  q+ 

0.0027 

0.0027 

0.0032 

2 

Maximum  Gr/Re. 

3.90  x I0“5 

4.84  x I0“5 

3.22  x I0"3 

Maximum  Mach  Number 

0.26 

0.25 

0.33 

x/D  for  Local  Bulk  Nusselt 
Numbers 

2.  1 -32.1 

2.1  - 32.1 

2. 1 -32.1 

with  axial  position  and  temperature  level  and  decreases  with  Reynolds  number. 

For  the  helium-argon  data  the  worst  situation  in  the  present  data  occurred 
~ 4 

at  Re^  ~ 3.1  x 10  and  M = 29.7  where  S^gg/^  reached  0.15  at  the  highest 
heating  rate;  thus,  for  this  condition  a ten  percent  uncertainty  in  heat 
would  cause  less  than  two  percent  uncertainty  in  the  heat  flux  to  the  gas. 

For  mixtures  with  lower  molecular  weights,  improved  convective  heat  trans- 
fer led  to  lower  heat  loss  ratios,  e.g.,  (q£  /q£)  ~ °-07  with  M = 15.3 

at  the  same  conditions. 

Adiabatic  friction  factors  were  measured  before  each  series  of 
heated  runs.  These  results  were  used  as  a check  of  the  measurements  of 
pressure,  mixture  molecular  weight,  and. flow  rate.  The  measured  friction 
factors  were  compared  to  the  experimental  correlation  of  Drew,  Koo,  and 
McAdams  Cl932U, 


f „ = 0.0014  + 0.125  Re 
DKM 


Figure  3 shows  the  measured  friction  factors  plotted  as  a function  of 
Reynolds  number.  Air  and  helium  data  points  are  included  for  comparison. 
All  the  measured  values  are  within  four  percent  of  correlation  (9) , and 
three-quarters  are  within  two  percent.  The  data  for  the  mixture  with 
M = 15.83  and  for  helium  show  the  best  agreement.  No  systematic  variation 
with  molecular  weight  seems  evident. 


2.5  Heat  Transfer  With  Constant  Properties 


In  order  to  deduce  the  turbulent  Prandtl  number  without  complica- 
tions introduced  by  its  possible  variation  with  heating  rate,  the  measure- 
ments were  extrapolated  to  the  constant  property  idealization  by  the  ap- 
proach of  Malina  and  Sparrow  Cl964].  For  this  method,  a series  of  experi- 
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and  McAdams  Correlation  for  Air,  Helium  and  Helium-Argon 
Mixtures. 


mental  runs  are  taken  with  the  same  inlet  Reynolds  number  and  gas  composi- 
tion but  with  successively  higher  heating  rates.  At  each  thermocouple 
location  the  measured  Nusselt  numbers  are  plotted  versus  the  local  tempera- 
ture difference,  T - T,  . An  extrapolation  to  T - T,  = 0 yields  the 
w b w b 

deduced  Nusselt  number  for  constant  property  conditions , Nu  . A compar- 

cp 

able  extrapolation  of  the  experimental  uncertainty  provides  an  estimate  of 
the  validity  of  these  data.  For  the  range  of  data  reported  the  uncer- 
tainty in  constant  properties  Nusselt  number  varied  from  nine  percent 
at  small  axial  distances  to  about  five  percent  at  large  distances. 

For  fully  developed  conditions  with  air  or  helium  as  the  fluid, 

the  ratio  Nu  /Nu  varied  from  0.94  to  1.0  without  any  evident  dependence 
Cp  DB 

on  Reynolds  number.  With  helium-argon  mixtures  of  Prandtl  number  0.42  to 
0.49  the  ratio  varied  from  0.83  to  0.93.  This  reduction  in  the  normalized 
value  corresponds  to  the  trend  predicted  in  Figure  2.  Since  the  reduction 
is  greater  than  predicted  by  the  analysis,  it  is  a first  indication  that 
a turbulent  Prandtl  number  greater  than  unity  may  be  appropriate  for  the 
mixtures . 

In  addition  to  equations  (1)  and  (2) , several  other  correlations 
have  been  recommended  for  heat  transfer  to  gases  with  fully  established 
temperature  profiles  downstream  in  tubes . Kays  C1966]  recommends  the 
relation 

Nu  = 0.022  Re°’8Pr°'6  (10) 

for  the  range  0.5  < Pr  < 1.0  for  the  thermal  boundary  condition  of  a con- 
stant wall  heat  flux.  This  analytical  correlation  agreed  with  the  present 
mixture  data  within  six  percent  so  it  can  be  considered  valid  to  Pr  = 0.42; 
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at  Pr  - 0.7  it  shows  slightly  better  agreement  than  the  Dittus-Boelter 
correlation  (2) 

Sleicher  and  Rouse  [1975]  suggested  a correlation  for  the  ranges 
0.1  < Pr  < 105  and  104  < Re  < 106, 


Ni^  = 5 + 0.015  Re™  Pr^  (11) 

where 

n = (1/3)  + 0.5  exp  (-0.6  Pr  ) 

w 


and 

m = 0.88  - 0. 24/(4  + Pr  ) 

w 

The  subscripts,  b,  f and  w,  refer  to  evaluation  of  properties  at  bulk,  film 
and  wall  temperatures,  respectively,  when  property  variation  is  significant. 

For  constant  properties  in  the  range  of  the  present  data  this  relation  pre- 
dicts lower  values  than  the  simpler  correlation  of  Kays.  At  Pr  ~ 0.7  the 

predictions  of  Sleicher  and  Rouse  correspond  to  the  lower  limit  of  the  data  and 

at  Pr  - 0.45  they  are  a few  percent  below  the  lower  limit  of  these  data. 

While  the  magnitudes  differ  from  the  data  slightly,  their  trend  with  Prandtl 

number  corresponds  to  the  present  predictions  and  measurement. 

In  order  to  deduce  the  turbulent  Prandtl  number  the  technique  of 
McEligot,  Picket  and  Taylor  [1976]  is  applied.  Essentially,  hypothesized 
distributions  Prfc(y)  are  introduced  in  the  numerical  analysis  and  the  pre- 
dicted Nusselt  numbers,  Nu(x),  are  compared  to  the  data,  Nu  (x) . Varia- 

cp 

tions  in  the  shape  of  Prt(y)  are  reflected  in  variation  in  the  shape  of 
Nu(x)  in  the  thermal  entry  region.  One  of  the  main  advantages  of  this  match- 
ing technique  is  that  it  concentrates  on  determining  Prfc  accurately  in  the 


region  which  is  most  important  for  predicting  surface  temperatures  in  flows 
heated  from  the  wall. 

In  this  study,  as  in  the  demonstration  by  McEligot,  Pickett  and 

Taylor  [1976],  the  experimental  uncertainty  at  small  axial  distances  was 

too  large  to  determine  d ?rt/d(y/rw)  clearly.  Consequently,  a constant 

value  of  Prfc  was  used  in  the  present  calculations;  their  sensitivity  tests 

showed  that  it  is  essentially  equivalent  to  an  effective  Pr  , for  the  wall 

region.  Only  results  for  x/D  greater  than  eight  were  used  to  determine  Prt  w- 

Figure  4 illustrates  examples  of  the  comparisons  between  measured 

Nu  (x)  and  predicted  Nu(x)  used  to  determine  Pr  . Examples  for  three 
op  t,w 

Reynolds  numbers  and  two  Prandtl  numbers  acre  shown.  Brackets  indicating 

the  estimated  experimental  uncertainties  of  the  Nu_^x)  are  also  included. 

Figure  4a  shows  the  measured  and  calculated  Nu  for  a helium-argon 

mixture  with  a molecular  weight  of  15.3,  Prandtl  number  of  0.42,  and 

Reynolds  number  of  32,000.  The  turbulent  Prandtl  number  is  seen  to  be 

about  1.1  or  slightly  less.  From  comparisons  on  similar  graphs,  Pr^ 

t,w 

was  estimated  to  be  1.1  0.1  for  helium-argon  mixtures  with  molecular 

weights  of  approximately  15,  Prandtl  numbers  of  0.42,  and  Reynolds  numbers 

between  32,000  and  55,200.  Measurements  and  predictions  are  shown  in 

Figures  4b  and  4c  for  a helium-argon  mixture  at  a molecular  weight  of  29, 

Prandtl  number  of  0.486,  and  Reynolds  numbers  of  31,600  and  82,100.  From 

similar  graphs,  Pr  was  estimated  to  be  1.0  + 0.1  for  the  range  of  helium- 
t,w  — 

argon  mixtures  with  molecular  weights  between  27  and  30,  Prandtl  numbers 
between  0.46  and  0.49,  and  Reynolds  numbers  between  31,600  and  102,000. 

The  effect  of  Reynolds  number  on  Pr.  ^ can  be  examined  qualitatively 
using  the  results  in  Figure  4b  and  4c,  but  the  experimental  uncertainty 
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A.  HELIUM -ARGON,  M = 15.3 


NUMERICAL  JeM  FROM  VAN  DRIEST/ REICHARDT 
ANALYSIS  1 CONSTANT  PROPERTIES 


Re  = 32000 
Pr  = 0.419 


B.  HELIUM-ARGON,  M =29.7 


Re=3lj600 
Pr  = 0.486 


Prt=  1. 1 - 
Prt*  1.2- 


C.  HELIUM  - ARGON,  M = 29.7 


Re  = 82J00 
Pr  =0.486 


. Prt=0.9 
;Prt=I.O 


,Prt  = 0.9 
,Prt=I.O 
Prt=  I . I 
. Prt=l.2 


Fiaure  4. 


Examples  of  the  Comparisons  Between  Measured  Nu  (x)  and  Predicted 
N u ( x ) as  Used  to  Deduce  Pr.,  for  Constant  Properties. 


prohibits  a quantitative  determination.  These  results  are  for  the  same 

Prandtl  number  (Pr  = 0.486) , but  two  different  Reynolds  numbers:  31,600 

and  82,100.  For  x/D  greater  than  eight,  and  at  the  lower  Reynolds  number, 

the  measured  Nu  (x)  data  are  slightly  below  the  calculated  Nu(x)  curve 
cp 

for  a Pr  of  1.0.  At  the  higher  Reynolds  number  and  same  axial  locations, 

1 1 w 

the  measured  Nu  (x)  are  slightly  above  the  prediction  for  Pr  =1.0. 

cp  t,w 

For  the  stated  conditions,  it  appears  that  Pr  may  have  a weak  dependence 

t f w 

on  Reynolds  number,  and  may  decrease  slightly  as  the  Reynolds  number  increases. 

The  effect  of  molecular  Prandtl  number  on  turbulent  Prandtl  number 

can  be  examined  using  the  observations  from  Figures  4a  and  4c  which  are 

summarized  in  Table  2.  The  air  data  of  McEligot,  Pickett  and  Taylor  Cl976] 

cam  also  be  considered  since  the  present  study  shows  Pr  to  vary  only 

t,w 

slightly  with  Reynolds  number.  For  the  range  shown  in  Table  2,  0.42  < Pr 
< 0.7,  the  turbulent  Prandtl  number  increases  as  the  molecular  Prandtl 
number  decreases  as  suggested  by  the  review  of  A.  J.  Reynolds  [jl975j. 

However,  reference  to  Figure  4 shows  that  predictions  based  on  O.  Reynolds' 
analogy  aire  still  valid  to  within  five  to  ten  percent  which  is  adequate 
for  mamy  engineering  applications. 


2.6  Heating  with  Property  Variation 


In  order  to  increase  power  densities  or  lower  the  weight  of  closed 
Brayton  cycle  systems* moderately  high  heating  rates  must  be  employed.  Con- 
sequently, the  temperature-dependence  of  the  fluid  transport  properties 
causes  3ignificiant  variation  in  properties  appearing  in  the  coefficients 


of  equations  (6)  and  coupling  of  the  equations.  The  results  and  correlations 
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Table  2. 

Variation  of  Pr,.  w 

t y W 

with  Respect  to  Molecular 

Prandtl  Number 

Gas 

Mol  ecu lar 
Weight 

Prandt! 

Number 

Pr 

t,w 

Reynolds 

Number 

Helium-argon  15.3 

0.419 

1 . 1 +0.1 

32,000 

Helium-argon  29.7 

0.486 

1 .0  + 0. 1 

31,600 

Air  * 

28.97 

0.72 

0.9  + 0. 1 

44,500 

*Mc£l igot. 

Picket  and  Taylor, 

[1976]  . 

i 

n 
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based  on  the  constant  property  idealizations  become  invalid.  In  this 
section  the  modifications  accounting  for  property  variation  are  examined 
and  the  application  of  the  turbulent  Prandtl  number,  determined  above  for 
constant  properties,  is  tested. 

Since  only  two  pressure  taps  were  attached  to  the  test  section, 
local  friction  factors  could  not  be  determined.  Overall  average  friction 
factors  with  heat  addition  were  compared  to  the  correlation  proposed  by 
Taylor  [1976] 

f = (0.0014  + 0.125  Re“°'32)  (T  /T.  )~°'5  (12) 

av  w w,av  b,av 


for  the  data  of  a wide  variety  of  experiments  with  gas  flow.  Most  previous 
measurements  agreed  within  ten  percent.  For  evaluation  of  this  expression, 
arithmetic  averages  of  the  conditions  at  the  two  pressure  taps  were  used  to 
calculate  the  average  temperatures  and  pressure  and,  hence,  density  and 
viscosity.  The  overall  friction  factor  was  determined  from  the  frictional 
pressure  drop. 


and  the  modified  wall  Reynolds  number  was  defined  as 


Re  = (GD.u  ) CT,  /T  ). 
w *w  b,av  w,av 


In  the  present  study  correlation  (12)  predicted  most  of  the  data 
within  four  percent;  all  are  within  ten  percent.  These  measurements  are 
presented  in  Figure  5.  Comparison  to  Figure  3 — which  uses  the  constant 
properties  limit  of  the  same  equation — shows  a general  increase  in  the 
normalized  data  with  heating.  This  increase  indicates  the  exponent  on  the 
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temperature  ratio  may  be  slightly  too  large,  leading  to  an  underprediction 


1 


of  the  average  wall  friction.  The  effect  appears  to  be  most  significant 
for  the  lower  molecular  weight  mixtures  in  low  Reynolds  number  flow. 

The  downstream  measurements  of  the  local  Nusselt  number  were  com- 
pared to  correlation  (11)  of  Sleicher  and  Rouse  [1975] . Though  intended  to 
account  for  property  variation,  the  values  predicted  were  15  to  40  percent 
lower  than  the  measurements  for  the  mixtures.  Thus,  this  correlation  could 
yield  excessively  conservative  estimates  of  surface  temperatures  in  design 
calculations. 

Mugee  [1968]  has  suggested  accounting  for  property  variation  and 
thermal  entry  development  by  applying  a multiplicative  factor  to  the  con- 
stant property  correlation.  In  the  range,  2.1  < x/D  < 81.6,  his  correlation 
predicted  97  percent  of  our  data  for  heat  transfer  to  air  and  helium  to  with- 
in 10  percent.  Since  the  transport  properties  of  helium-argon  mixtures  vary 
with  temperature  in  approximately  the  same  manner  as  for  air  and  helium,  the 
same  exponent  was  taken  for  the  temperature  ratio.  Kays  [1966]  discusses 
the  effect  of  Prandtl  number  variation  on  the  thermal  entry  behavior  for 
circular  tubes.  As  the  Prandtl  number  decreases  the  thermal  entry  region  is 
expected  to  become  more  pronounced,  so  a larger  coefficient  is  appropriate 
for  the  axial  development  term.  These  considerations  plus  a correlation  of 
the  constant  properties  results  [Pickett,  1976]  lead  to 

Nu,o  = 0.021  Re°‘8  Pr^‘  55  [ (T  y^)"0'4  + 0-85  D/x]  (13) 

which  predicts  92  percent  of  the  present  helium/argon  data  to  within  ten 
percent  for  the  range  2.1  < x/D  < 81.6.  The  greatest  discrepancy  occurs  in 
the  range  4 < x/D  < 16  at  high  heating  rates;  in  this  situation  the  measure- 
ments are  underpredicted  by  5 to  15  percent. 
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The  validity  of  the  deduced  values  of  the  turbulent  Prandtl  number 


were  tested  for  the  conditions  of  two  experimental  runs,  one  with  slight 
property  variation  and  the  other  with  strong  heating.  Both  runs  were  with 
the  same  mixture  with  M = 29.7  and  Pr  ~ 0.49  and  essentially  the  same  inlet 
Reynolds  number.  For  predictions  the  measured  wall  heat  flux  distribution 
was  used  for  the  thermal  boundary  condition  and  exponents  in  equations  (5) 
were  taken  as  a = 0.772  and  b = 0.741  to  account  for  property  variation. 

From  the  constant  property  results  the  turbulent  Prandtl  number  was  taken 
as  1.02. 

Comparisons  between  the  two  heating  rates  and  the  numerical  predic- 
tions are  presented  as  Figure  6.  These  demonstrate  that  the  combined  effect 
of  the  strong  heating  rate  and  the  property  variation  is  a reduction  in 
Nu^  of  about  1/3  in  the  downstream  region.  The  difference  in  Nu  decreases 
as  the  thermal  entry  is  approached  while  the  temperature  ratio  peaks  near 
ten  to  fifteen  diameters,  the  region  where  correlation  (13)  shows  the 
greatest  discrepancy. 

The  numerical  predictions  and  data  agree  well  for  both  heating 
rates.  This  observation  is  to  be  expected  for  the  lower  heating  rate, 
q"1"  = 0.0006,  since  the  maximum  difference  between  wall  and  bulk  properties 
was  about  eleven  percent.  Thus,  the  conditions  did  not  differ  substantially 
from  the  constant  properties  idealization  under  which  the  turbulent  Prandtl 
number  was  deduced. 

At  the  higher  heating  rate,  q+  = 0.0032,  the  property  variation  at 
a cross  section  approaches  sixty  percent.  Thus,  it  is  reasonable  to  ask 
whether  the  heating  causes  a significant  effect  on  the  turbulent  Prandtl 
number  in  the  important  wall  region.  Apparently,  it  does  not.  Applying 


Figure  6.  Comparison  of  data  to  Numerical  Predictions  Accounting  for  Transport  Property 
Variation. 


the  value  deduced  for  constant  properties  leads  to  close  agreement  between 
predicted  and  measured  Nusselt  numbers  it  this  heating  rate  as  well  (Figure 
6) . While  not  an  exhaustive  examination,  this  test  is  an  initial  indica- 
tion that  Pr  is  not  a strong  function  of  heating  rate.  Extension  to  other 
molecular  Prandtl  numbers,  Reynolds  numbers  and  heating  rates  is  the  subject 
of  further  work  now  in  progress. 

3.  NUMERICAL  PREDICTION  OF  WALL  FRICTION 
IN  LAMINAR  FLOW  WITH  PROPERTY  VARIATION 


3.1.  Background 

In  gas  turbine  cycles  the  regenerative  heat  exchanger  is  typically 
constructed  of  parallel  plates,  with  short  fins  attached  forming  additional 
parallel  surfaces.  Consideration  of  the  heat  transfer  performance  versus 
pumping  power  requirements  of  these  heat  exchangers  usually  results  in  design 
for  operation  in  the  laminar  or  transitional  flow  regime.  For  laminar  flow 
heat  exchangers , the  streamwise  length  of  the  fins  is  shortened  in  order  to 
take  advantage  of  the  increased  heat  transfer  coefficient  of  developing 
boundary  layers  by  continually  reinitiating  the  boundary  layer.  The  thermal 
boundary  condition  is  an  approximately  constant  wall  heat  flux.  As  a guide 
to  the  effects  of  mixture  composition  and  property  variation  in  such  geo- 
metries, the  present  work  - for  a first  objective  - investigates  the  simul- 
taneous development  of  laminar  thermal  and  velocity  boundary  layers  in  the 
entry  region  of  parallel  plate  ducts. 

The  previous  report  [Taylor  et  al. , 1976]  (a)  briefly  discussed 
related  work  which  can  be  extended  to  provide  improved  guidance  to  the 
designer  of  regenerative  heat  exchangers  for  mixtures  of  noble  gases, 


(b)  summarized  pertinent  knowledge  of  their  transport  properties  and 
demonstrated  the  generalizations  possible  to  reduce  the  analytical  task, 

(c)  outlined  the  numerical  solution  of  the  governing  equations,  and  (d) 
presented  the  heat  transfer  results  of  interest  in  design  in  terms  of 
lengthwise  mean  parameters.  The  numerical  results  are  summarized  in  the 
present  report  as  Table  3.  The  mean  Nusselt  number,  evaluated  in  terms 

of  properties  evaluated  at  the  average  bulk  temperature,  could  be  approxi- 
mated as 

Nu  = [8.2352  + 1.9312/(Pr°'254L*)]1/2  (14) 

cl 

This  equation  is  of  the  order  of  five  percent  lower  than  the  numerical 
prediction  near  L*  = 0.01.  For  moderate  heating  rates  it  would  be  within 
about  ten  percent  of  the  numerical  results  and  would  be  low;  in  heat  ex- 
changer design  this  would  lead  to  units  slightly  longer  than  necessary. 
Alternatively,  the  constants  in  equation  (14)  could  be  optimized  for  another 
range  at  the  expense  of  the  accuracy  of  predictions  in  the  immediate  entry. 

It  was  found  that  the  use  of  properties  evaluated  at  the  average  film  temp- 
eratures did  not  collapse  the  variable  properties  predictions  to  the  constant 
properties  result  as  well  as  the  bulk  temperature  did. 

3.2.  Prediction  of  Wall  Friction  wi th  Property  Variation 

While  most  analyses  presently  available  for  developing  flows  present 
friction  results  in  terms  of  the  wall  shear  stress  evaluated  from  the  velocity 
gradient  at  the  wall  (f  } , this  approach  is  not  of  use  to  the  designer  when 
the  velocity  profile  is  changing  substantially,  as  in  the  entry  or  when  a gas 
is  heated.  To  predict  the  required  pressure  drop  with  a one-dimensional 


Q*  - 0,  constant  properties 
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design  procedure,  one  uses  the  "apparent"  friction  factor,  f , 

ap 

the  wall  shear  determined  by  treating  the  momentum  change  as  one-dimensional. 

The  same  treatment  is  often  employed  in  experiments  where  size  prohibits 

velocity  profile  measurements.  Both  methods  of  presentation  can  be  chosen 

with  numerical  results;  consequently,  Bankston  and  McEligot  [1970]  were  able 

to  demonstrate  (a)  the  numerical  values  of  f and  f can  differ  substantially 

s ap 

and  (b)  discrepancies  earlier  thought  to  exist  between  experiments  and  an- 
alyses were  primarily  due  to  the  differences  in  the  definition  used  for  the 
friction  factors. 

As  with  the  heat  transfer  results  we  concentrate  in  presenting  a 
mean  apparent  friction  factor, 


f 

a 


^h 

4L 


£ 

G2/2gc 


A (p 

0 L 


Pbgc 


(15) 


(The  local  apparent  friction  factor,  f^,  appearing  in  Table  3 is  defined 

in  the  analogous  derivative  form  with  d/dx  replacing  (1/L)  A ).  When  , and 

0 L D 

constant,  the  second  term  in  brackets  does  not  change  and  the  definition 
reduces  to  that  of  Shah  and  London  [1971] . With  constant  fluid  properties 
one  solves  the  flow  problem  only,  so  the  result  is  independent  of  Prandtl 
number  and  can  be  written  as  a single  function  fa(L+)  which  approaches 
f * Re/24  as  L+  becomes  large.  This  function  may  be  found  tabulated  in 
Table  3 or  can  be  derived  from  earlier  local  results  [Schade  and  McEligot, 
1971;  Bodoia  and  Osterle,  1961] . For  a continuous  approximation,  the  approach 
of  Schlunder  [1974]  can  be  used  as  in  the  heat  transfer  results  to  give 

f • Re/24  = /I  + 0.0788/L+  (16) 
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which  represents  the  numerical  results  well  in  the  immediate  entry  but  is 
4 to  5 percent  high  in  the  range,  0.05  < L+  < 0.2. 

With  varying  transport  properties,  the  energy  equation  is  coupled 
to  the  flow  equations  via  the  temperature-dependent  viscosity  and  density, 
so  the  wall  friction  also  becomes  a function  of  the  Prandtl  number  and  the 
heating  rate.  Again  the  question  arises  as  to  the  better  method  of  account- 
ing for  the  fluid  property  variations.  Predictions  of  friction  are  not  as 
well  behaved  as  heat  transfer  parameters.  In  contrast  to  the  heat  transfer 
results,  direct  use  of  the  average  bulk  properties  in  f , Re  and  L+  does  not 

cl 

collapse  the  results  nicely  around  the  prediction  based  on  constant  properties; 

the  main  effect  is  to  spread  the  curves  towards  larger  L.  + as  Q+  increases. 

ba 

The  effect  of  heating  rate  on  apparent  wall  friction  is  presented 
in  Figure  7 partially  in  terms  of  average  bulk  properties.  That  is,  p^a  is 
used  for  the  coefficient  in  equation  (15)  and  Re^  is  defined  as  before  but 
the  non-dimensional  length  is  based  on  inlet  properties,  i.e.,  L*.  With  this 
representation,  heating  increases  f^.  Re^  considerably  more  than  Nu^a  is 
raised  at  the  same  level  of  Q*.  At  lengths  greater  than  L+  = 0.1  the  curves 
with  heating  approach  the  constant  properties  curve  only  slowly,  although 
T^/T  is  close  to  unity,  as  the  heated  entry  continues  to  affect  the  inte- 
grated results  far  downstream.  Close  inspection  of  the  trends  for  the  high- 
est heating  rates  shows  that  as  Pr  increases  a convergence  - from  heated 
entry  behavior  towards  agreement  with  constant  property  behavior  - is  moved 
further  downstream.  This  effect  corresponds  to  the  thermal  boundary  layer 
and  shear  boundary  layer  as  the  Prandtl  number  changes:  Nu(L*)  shows  only 
a moderate  effect  of  Pr,  so  for  the  same  heating  rate  the  value  of  T^/T^  is 
almost  the  same  at  equal  values  of  L*  rather  than  L+,  thus  T^/T^  approaches 
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unity  for  Pr  = 0.2  at  earlier  values  of  L than  for  Pr  = 2/3  and  the 

variation  of  properties  across  the  channel  is  less  for  Pr  = 0.2  at  the 

same  L+.  While  the  friction  predictions  for  heating  approach  the  adiabatic 

prediction  as  increases,  those  for  cooling  diverge;  this  result  also 

corresponds  to  the  trend  of  property  variation  since  the  ratio  T,  /T 

ba  wa 

increases  downstream  for  cooling  as  described  earlier  in  the  section  on 
heat  transfer. 

As  with  the  heat  transfer  results,  the  apparent  effect  of  property 
variation  on  wall  friction  is  sensitive  to  the  choice  of  reference  tempera- 
ture. With  average  film  temperature  for  the  reference,  the  shape  of  the 
resulting  curves  differs  from  the  shape  with  bulk  temperature  as  reference. 
In  Figure  8 the  product  ffa  • Refa/24  is  plotted  against  L*;  pfa  is  used 

in  the  coefficient  in  equation  (15)  and  Re„  is  based  on  . There  is  no 

fa  fa 

advantage  in  comparison  on  the  basis  of  L*a  since  results  are  shifted  then 
further  to  the  right  (with  heating)  so  that  for  L+  > 0.01  the  difference 
from  the  adiabatic  prediction  is  increased.  For  heating:  the  friction 
parameters  are  reduced  for  short  lengths;  then  the  predictions  converge  with 
and  across  the  constant  properties  curve  and  remain  slightly  greater  at 
larger  distances.  In  comparison  to  the  bulk  property  predictions,  the  ef- 
fects for  short  and  long  ducts  are  approximately  the  same  magnitude  with 
strong  heating,  but  for  intermediate  lengths  and  for  Q+  < 10  a display  in 
terms  of  film  properties  shows  signicantly  less  variation.  In  the  range 
-2  < Q+  < 2 there  is  no  significant  effect  of  heating  until  approaches 
0.1  with  film  properties  and  then  the  effect  is  only  of  the  order  of  five 
to  ten  percent.  As  is  the  case  for  bulk  properties,  the  convergence  towards 
the  adiabatic  prediction  is  at  successively  greater  distances  (L*)  as  the 


Prandtl  number  increases , but  for  the  same  condition  it  is  several  times 
earlier  with  film  properties.  With  cooling:  the  directions  of  the  trends 
are  reversed  but  for  Q+  > -2  they  are  essentially  again  negligible  for 
entry  problems. 

It  is  not  clear  from  Figure  7 and  B which  approach  is  better: 
property  ratio  or  film  temperature.  The  property  ratio  approach  would  be 
represented  as 


(f. 
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Re  (L  ))/(f 

O dL 


Re  (L+)  ) = (T  /T.  )' 

o cp  wa  ba 
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so  this  quotient  is  plotted  versus  temperature  ratio  in  Figure  9 to  examine 

the  suitability  of  a single  exponent.  A complicated  pattern  appears.  In 

contrast  to  the  expectation  of  Kays  and  London  [1964] , the  general  trend  is 

a substantial  increase  with  temperature  ratio.  There  are  slight  differences 

with  Prandtl  number  but  the  trends  are  mostly  the  same.  An  exponent  q of 

the  order  of  unity  would  overpredict  the  friction  factor  at  the  higher 

temperature  ratios  and  underpredict  it  at  lower  values.  For  cooling,  q = 1 

is  valid  within  a few  percent.  For  moderate  heating,  the  necessary  value 

of  q (i.e.,  the  slope  of  a line  from  the  origin  on  this  logarithmic  plot) 

varies  with  length  L*:  it  is  approximately  constant  as  the  temperature 

ratio  increases  with  length  then  increases  gradually  as  the  ratio  drops 

for  successively  longer  ducts.  The  latter  effect  is  a consequence  of  the 

slow  convergence  of  f • Re^  to  the  adiabatic  curve  for  long  ducts  as 

discussed  earlier.  It  is  seen  that  a function  q(L+,  Q+,  Pr)  would  be 

o 

necessary  to  describe  the  detailed  behavior.  For  Q+  = 2 and  L*  < 0.6,  an 
exponent  q = 1.5  would  reduce  the  difference  from  the  constant  properties 
curve  from  13  percent  to  a 7 percent  discrepancy.  With  Q+  = 10,  q = 1.2 
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Figure  9.  Examination  of  Property  Ratio  Method  for  Mean  Wall  Friction. 


is  a better  approximation,  but  the  discrepancy  would  still  reach  twenty 
percent.  These  comments  and  comparison  of  Figures  7 and  9 suggest  that 
the  two  methods  have  approximately  the  same  overall  accuracy  for  Q+  < 2 
with  a slight  advantage  to  the  film  properties  approach  for  short  ducts. 

For  moderate  heating  - to  Q = 10  - the  film  property  method  is  clearly 
superior,  while  at  higher  heating  rates  both  methods  show  regions  where 
the  simple  correlations  would  mislead  the  designer  substantially. 

It  is  perhaps  inconvenient  for  the  designer  to  have  one  method 
perform  better  for  heat  transfer  while  the  other  is  preferable  for  wall 
friction,  but  the  difficulty  should  be  negligible  provided  the  present 
definitions  of  the  parameters  are  used.  Once  the  heat  transfer  problem 
is  solved  for  the  wall  temperature  using  average  bulk  properties,  the 
average  film  temperature  can  be  calculated  from  the  results  and  can  then 
be  employed  to  predict  the  wall  friction  behavior. 

Analytical  correlations  such  as  equations  (14)  and  (16)  are  useful 
for  parameter  studies  of  systems  and  for  initial  sizing  of  components  when 
hundreds  to  thousands  of  individual  configurations  may  be  calculated.  When 
greater  accuracy  is  needed  in  final  design  decisions  - or  if  variable  wall 
heat  flux  should  be  treated  - the  numerical  analysis  can  be  employed  directly. 
With  the  direct  application  of  the  program,  the  question  of  definitions  of 
the  non-dimensional  parameters  is  avoided;  the  engineer  can  choose  definitions 


to  suit  his  own  convenience,  including  direct  presentation  in  temperatures, 
pressure  and  lengths  in  units  of  his  choice. 


4.  INITIAL  MEASUREMENTS  WITH  CLOSED  LOOP 


The  schematic  description  of  the  closed  loop  for  measurements  with 
expensive  gas  mixtures,  such  as  helium  and  xenon,  has  been  presented  by 
Taylor  et  al.  [1976] . The  basic  loop  apparatus  was  completed  during  the 
summer,  1976.  Thermal  instrumentation  has  been  interfaced  to  our  PDP-8M/ 
HP3480-5  digital  data  acquisition  system  [McEligot,  1975b] . All  components 
of  this  system  have  been  received  except  the  KL8J-A  interface  for  direct 
connection  to  the  University  DEC-lO/CDC  6400  system;  this  interface  is  on 
order.  In  the  present  configuration  data  stored  in  the  PDP-8M  can  be  trans- 
ferred on  DEC  tape  to  the  DEC-10  by  courier.  The  initial  data  acquisition 
programs  are  in  operation. 

The  first  test  section  in  the  loop  was  a circular  tube  of  0.253 

«• 

inch  0D  and  0.0115  inch  thick  walls,  with  60  diameters  between  electrodes. 
Heat  loss  and  resistivity  calibrations  were  conducted  from  room  temperature 
to  over  1200°F.  Adiabatic  friction  measurements  spanned  the  range  35,000 
< Re  < 83,000  with  air  in  order  to  test  flow  measurement  techniques  pri- 
marily. Preliminary  tests  included  the  following  ranges: 

Molecular  Prandtl  Reynolds  Max-Wall 


Gci  s 

Weiqht 

Number 

Number 

Temperature 

System 

Argon 

40.0 

0.67 

7,600 

510°F 

rw/Diaphr 2 
Closed  Loop!  /pump 

Air 

29.0 

0.72 

30,000-80, 

000  1010°F 

Open  Loop 

He-Ar 

31.0 

0.50 

40,000 

230°F 

. _ ,w/Piston 

Closed  Loop!  piJmp 

Air 

29.0 

0.72 

50,000 

200°F 

Closed  Loop{w/,p;'-ston 
Pump 

Mixtures 

of  gases 

other  than 

the  noble  gases 

are  expected  to  pro- 

vide  reduced  molecular  Prandtl  numbers  as  well.  Since  many  common  gases 
are  less  expensive  than  the  noble  gases  their  mixtures  are  being  investi- 
gated. The  properties  of  hydrogen-carbon  dioxide  mixtures  are  currently 
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being  surveyed.  The  estimated  variation  in  Prandtl  number  is  plotted  on 
Figure  10  for  this  mixture,  it  is  seen  that  the  minimum  is  Pr  ~ 0.34  at 
M - 15.  A supply  at  this  molecular  weight  and  at  M ~ 29  has  been  obtained 
for  comparison  to  the  data  of  section  2 and  to  extend  our  data  to  lower  Pr. 

On  March  3,  during  a preliminary  test  with  helium  flow  at  a moder- 
ate heating  rate  the  test  section  failed  in  the  vicinity  of  the  pressure 
tap  attachment  at  x/D  = 56.  Inspection  showed  evidence  of  brittle  fracture 
along  grain  boundaries  as  well  as  some  molten  regions  It  is  hypothesized 
that  the  tube  failed  due  to  intergranular  fracture . and  that  as  the  two 
sections  parted  the  electric  arc  caused  local  melting.  Failure  analysis 
with  scanning  electron  microscope  and  microprobe  of  the  University's  Space 
Sciences  Center  is  planned.  A new  test  section  is  under  construction  and 
should  be  installed  by  mid  April. 

5.  CONVECTION  HEAT  TRANSFER  FROM  ROUGH  SURFACES 


Proposals  for  increasing  the  performance  of  current  Brayton  systems 
include  operation  at  higher  pressure  levels,  use  of  alternate  fluids, 
increased  use  of  roughness  elements  and  extended  surfaces,  and  selection 
of  materials  capable  of  higher  temperatures.  The  highest  temperatures 
are  in  the  heating  component.  The  proposed  Gas  Cooled  Fast  Breeder  Reactor 
(GCFBR)  and  the  defunct  nuclear  rocket  propulsion  program  faced  comparable 
requirements  for  high  power  densities.  The  GCFBR  technology  is  planning 
on  using  helium  at  high  pressure  with  artificially  roughened  coolant  chan- 
nels fabricated  from  high  temperature  alloys.  On  the  other  hand,  the  Nerva 
program  chose  graphite,  pyrolitic  carbon  and  niobium  ca.rnide  to  operate 
at  higher  temperatures  than  possible  with  metal  alloys.  Other  proposals 
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suggest  using  ceramics  such  as  silicon  nitride  (Si3N4)  for  the  heater 
and  turbine  blades.  Each  of  these  approaches  potentially  involves  basic 
convection  problems  which  have  not  been  significant  in  existing  small 
Brayton  engines  for  Naval  applications. 

Our  studies  to  date  have  concentrated  on  examining  convective  heat 
transfer  and  wall  friction  for  alternate  working  fluids,  the  noble  gases 
and  their  mixtures.  The  expense  of  the  high  molecular  weight  Noble  gas 
mixtures  leads  to  concern  over  the  possibility  of  leakage  in  large  scale 
plants.  While  such  leakage  is  likely  to  be  primarily  the  lower  molecular 
weight-  and  cheaper-component,  other  mixtures  of  more  common  gases  could 
be  used  to  provide  some  of  the  advantages  of  the  Noble  gas  mixtures  at 
lower  expense.  Accordingly,  we  have  begun  measurements  with  alternate  gas 
mixtures  as  described  in  the  previous  section. 

Development  of  ceramic  tubes  for  the  heating  component  is  critical 
to  the  success  of  the  some  proposed  improved  Brayton  cycles.  Likewise, 
ceramic  turbine  blades  are  suggested  as  a means  of  increasing  turbine  in- 
let temperatures  without  the  need  for  internal  blade  cooling  [Johansen  and 
Wallace,  1976] . Fabrication  studies  show  a rough  surface  with  closely- 
spaced,  approximately  triangular,  ribbed  shape;  the  strength  of  blade  de- 
pends strongly  on  the  orientation  of  these  ribs.  With  heating  in  an  oxi- 
dizing atmosphere,  the  surface  texture  changes  considerably  in  a few  hours; 
one  would  expect  a comparable  variation  during  the  operating  life  of  an 
inert  gas  system  at  a low  concentration  of  impurities  such  as  residual  air 
from  initial  changing  of  the  system.  The  change  in  texture  will  modify 
the  frictional  pressure  drop  and  convective  heat  transfer  characteristics 
of  the  ceramic  surface. 
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The  estimated  geometry  of  some  typical  roughened  surfaces  is  pre- 
sented in  Table  4.  Currently  under  design  is  a test  section  for  measuring 
velocity  and  temperature  distributions,  heat  transfer  parameters  and  pres- 
sure drop  for  surfaces — representing  both  naturally  occurring  ceramic  sur- 
faces and  artificially  roughened  superalloys — at  heating  rates  causing 
significant  property  variations  as  for  high  power  densities. 

To  increase  power  densities  in  the  heater  of  a Bray ton  system  for 
lightweight  ship  propulsion,  the  entire  length  of  the  coolant  channel  does 
not  need  to  be  roughened.  For  a thermal  boundary  condition  of  a specified 
heating  rate,  q^(x) , the  tube  wall  temperature  peaks  in  the  downstream 
region  of  the  channel.  By  using  a rough  surface  only  in  the  section  with 
maximum  wall  temperatures,  the  peak  value  can  be  reduced  below  limits  im- 
posed by  material  limitations  and/or  the  heating  rate  can  be  raised  so  the 
turbine  inlet  temperature  is  increased.  The  pressure  drop  is  increased 
only  in  the  vicinity  of  the  rough  surface  so  substantially  larger  heat 
transfer  improvement  can  be  accommodated  locally  without  roughening  the 
entire  channel  (and,  thus,  increasing  the  pressure  drop  throughout). 

Basic  data  necessary  for  design  include  the  development  of  the 
internal  thermal  and  momentum  boundary  layers  caused  by  the  change  in  sur- 
face texture,  the  fully  developed  Stanton  number  and  friction  factor,  and 
the  redevelopment  to  a smooth  wall  boundary  layer,  all  with  heating  rates 
causing  significant  property  variation.  Currently  data  are  required  for  a 
variety  of  shapes,  but  advanced  numerical  methods  for  flows  with  recircula- 
tion may  provide  means  to  reduce  experimentation  in  the  future. 

The  experimental  apparatus  (Figure  11)  is  to  use  replaceable  test 
sections  so  that  various  forms  of  roughness  elements  can  be  tested.  The 
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Table  4.  Characterization  of  Rough  Surfaces  (Typical  Design  Geometry  - All  Dimensions  in  Inches). 


Figure  1 1 . --continued 


test  section  will  consist  of  two  electrically  heated  parallel  plates  about 

twelve  inches  wide,  spaced  one  inch  apart.  Parallel  plates  are  chosen  to 

allow  use  of  the  laser  Doppler  anemometer  for  velocity  measurement;  one 

inch  plate  spacing  is  chosen  to  reduce  natural  convection  effects  in  the 

initial  measurements.  A sketch  of  the  initial  heated  plate  is  shown  in 

Figure  12.  Initial  elements  will  be  square  ribs  to  extend  existing  data 

to  examine  the  effect  of  .heating  rate  through  the  temperature  dependence 

of  the  properties.  Pitch-to-height  ratio,  p/h,  will  be  of  the  order  of 
. + 

ten  with  h about  20.  A schematic  diagram  of  the  test  section  instrumen- 
tation has  been  shown  earlier  [McEligot,  1975a]  . Air  will  be  used  as  the 
gas  (Pr  - 0.7).  Wall  thermocouples  and  the  power  measurements  provide  the 
information  to  calculate  local  Stanton  numbers.  Temperature  profile  measure- 
ments will  be  made  with  a commercial  hot  w-ire  anemometer  operated  as  a re- 
sistance thermometer  to  determine  thermal  boundary  layer  development  and 
the  heat  transfer  roughness  function, 

G = t+  * *n(y/h) 

The  laser  Doppler  anemometer  will  be  used  to  measure  the  axial  velocity 
profile  to  deduce  the  momentum  boundary  layer  development  and  the  momentum 
transfer  roughness  function, 

R = u+  - ^£n(y/h) 

The  roughness  functions  provide  the  information  necessary  to  set  the 
boundary  conditions  to  solve  the  energy  equation  and  momentum  equation 
numerically  for  design  predictions. 
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Figure  12.  Initial  Test  Surface  (to  scale! 
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Currently,  negotiations  with  fabricators  are  in  progress  for  con- 
struction of  the  test  surface  and  materials  are  on  order.  The  single- 
component, laser  Doppler  anemometer  system  is  complete  with  the  exception 
of  the  TSI  1090  frequency  tracker  which  is  expected  to  be  delivered  in 
April.  The  receiving  optics  from  OEI  (Karlsruhe!  were  tested  by  the 
Principal  Investigator  at  the  Universitat  Karlsruhe  during  July  1976 
prior  to  shipment. 

6.  CONCLUSIONS 

The  general  goal  of  the  project  is  to  develop  the  ability  to 
predict  accurately  temperature  and  velocity  distributions,  pressure  drop 
and  wall  heat  fluxes  in  components  of  the  closed  Brayton  cycle.  Since 
mixtures  of  inert  gases  offer  advantages  when  used  as  the  working  fluid, 
the  main  emphasis  of  the  present  study  has  been  to  investigate  the  ques- 
tion whether  existing  design  correlations  for  pure  gas  flows  can  be  used 
reliably  for  such  mixtures. 

Examination  of  the  properties  of  inert  gas  mixtures  shows  that, 
for  concentrations  providing  improved  heat  transfer  performance,  the 
Prandtl  number  is  reduced  compared  to  pure  gases.  Heat  transfer  data 
have  not  previously  been  available  in  this  range  (Pr  < 0.5)  so  the  effect 
of  the  Prandtl  number  variation  had  not  been  measured.  Since  high  power 
densities  require  high  wall  heat  fluxes,  large  temperature  differences 
may  be  encountered:  calculations  of  mixture  thermal  conductivity  and  vis- 
cosity vary  approximately  as  the  0.7  to  0.8  power  of  the  absolute  tempera- 
ture . 
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Operating  conditions  for  heat  transfer  components  of  Bray ton  cycle 


systems  for  Naval  applications  are  expected  to  be: 

Regenerative  heat  exchangers  plate/fin  geometry,  laminar  or 

transitional  flow,  approximately 
constant  wall  heat  flux 

Heater  tubes  or  ducts,  turbulent  flow, 

specified  wall  heat  flux  distri- 
bution 

Cooler  plate/fin  geometry,  laminar  or 

transitional  flow,  approximately 
constant  wall  temperature 

To  simulate  conditions  in  the  heater,  numerical  analysis  was 
applied  to  turbulent  flow  in  a circular  tube  with  constant  wall  heat  flux. 
For  constant  fluid  properties,  it  was  predicted  that  the  local  Nusselt 
number  would  be  reduced  below  the  value  given  by  the  existing  Dittus- 
Boelter  correlation  and  Colburn  analogy  for  gases;  the  reduction  would 
become  greater  as  the  Prandtl  number  is  decreased  and  Reynolds  number  is 
increased.  These  trends  were  confirmed  by  the  experiments  described  below. 

While  numerical  predictions  for  laminar  flow  require  no  empirical 
assumptions,  turbulent  flow  predictions  are  based  on  the  assumption  of  a 
turbulent  Prandtl  number.  However,  its  dependence  on  molecular  Prandtl 
number  is  not  yet  known  so  the  turbulent  predictions  must  be  tested  by  ex- 
periment. Measurements  of  wall  temperature  and  pressure  drop  were  obtained 
in  a vertical,  resistively  heated,  circular  tube  with  a nominal  diameter  of 
1/8  inch  (3mm)  with  a 92  diameter  adiabatic  entry  followed  by  a 98  diameter 
heated  section.  Gases  were  air,  helium  and  helium-argon  mixtures  of  mole- 
cular weights  approximately  15.8  and  29.7.  The  range  of  variables  covered 
was  Reynolds  numbers  from  30,000  to  102,000,  Prandtl  number  from  0.42  to 
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0.7  and  wall  temperatures  from  room  temperature  to  1500°R,  wall-to-bulk 

temperature  ratios  up  to  1.9  and  axial  distances,  2.1  < x/D  < 32. 

The  data  for  heat  transfer  parameters  were  extrapolated  to  a 

wall-to-bulk  temperature  difference  of  zero  to  approximate  the  constant 

properties  idealization.  Comparison  between  these  deduced  data  and 

numerical  predictions  based  on  hypothesized  values  of  Pr  , the  turbuler.* 

tw 

Prandtl  number  near  the  heated  wall,  provided  estimates  of  the  values  of 

Earlier  tests  had  shown  that  for  surface  heating  the  resulting  wall 

temperatures  are  primarily  sensitive  to  Pr  and  that  values  of  Pr  in  the 

tw  t 

turbulent  core  are  not  particularly  important.  The  turbulent  Prandtl 
number  in  the  wall  region  was  deduced  to  be  1.1  +_  0.1  for  Pr  ~ 0.42, 

1.0  +.0.1  for  Pr  - 0.49  and  0.9  +.0.1  for  Pr  ~ 0.7  for  constant  properties; 
variation  with  Reynolds  number  was  slight.  These  results  do  not  differ 
substantially  from  Reynolds  analogy  which  is  commonly  employed  in  analyses. 

With  heating  at  conditions  which  caused  property  variation  to  be 
significant,  the  helium-argon  data  were  correlated  as 

0.3  0.55  -0.4 

Nu  = 0.021  Re  Pr  [ (T  /T,  ) + 0.85  D/x] 

b b w b 

for  0.42  < Pr  < 0.5  to  within  ten  percent  for  92  percent  of  the  data.  The 
overall  pressure  drop  data  were  correlated  in  terms  of  the  average  friction 
factor 

f = (0.0014  + 0.125  Re"0'32) (T  /T  ) “° ‘ 5 
av  w w,av  b,av 

The  numerical  analysis  was  extended  to  account  for  fluid  property 

variation.  For  predictive  purposes  the  value  of  Pr  was  taken  from  the 

tw 
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results  of  the  constant  property  study.  Comparison  to  measurements  with 

a helium-argon  mixture  of  Pr  = 0.49  showed  the  predictions  to  be  adequate 

at  (T  /T  ) ~ 1.8  and  Re.  ~ 3 x 10. 

w b max  l 

To  simulate  conditions  in  the  regenerative  heat  exchanger,  numeri- 
cal analysis  was  applied  to  laminar  flow  between  parallel  plates  with  con- 
stant wall  heat  flux.  Our  heat  transfer  results  were  reported  earlier. 

For  low  temperature  differences  such  that  the  fluid  properties  are  approxi- 
mately constant,  the  existing  Sieder-Tate  correlation  was  found  to  under- 
predict the  Nusselt  number  substantially,  to  represent  the  Prandtl  number 
dependence  incorrectly  and  to  vary  with  non-dimensional  length  L*  in  a 
different  manner  than  the  numerical  predictions.  A new  correlation  was 
obtained  for  the  mean  Nusselt  number  for  simultaneous  velocity  and  thermal 
boundary  layer  development, 

Nu  = [3.2352  + 1.9312/(Pr°'254  L*)]i/2 

cl 

This  correlation  agrees  with  the  predictions  to  within  about  five  percent. 
In  the  entry  region  the  mean  Nusselt  number  increases  about  17%  as  the 
Prandtl  number  is  reduced  to  0.2  from  2/3.  It  was  found  that  when  strong 
heating  causes  the  fluid  properties  to  vary  significantly,  the  mean  Nusselt 
number  is  only  increased  slightly  above  the  value  predicted  by  the  constant 
property  analysis  provided  the  properties  are  evaluated  at  the  average  bulk 
temperature.  If  properties  are  evaluated  at  the  average  film  temperature, 
the  change  in  Nusselt  number  is  greater  and  appears  as  a reduction  rather 
than  an  increase. 

For  wall  friction  under  the  constant  properties  idealization,  the 
average  friction  factor  predicted  by  the  numerical  analysis  can  be  approxi- 
mated by 
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+ 1/2 

(f  • Re/24)  = [1  + 0.0788/L  ] ' 

d 


to  within  ten  percent  for  L > 0.001.  In  contrast  to  the  heat  transfer 
results,  the  bulk  properties/property  ratio  method  does  not  collapse  the 
friction  predictions  well.  However,  in  the  range  -2  < Q+  < 10  the  film 
temperature  approach  for  evaluating  fluid  properties  agrees  with  the  con 
stant  properties  friction  relation  above  so  it  is  recommended. 
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